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Abstract. We study scalar field and electromagnetic perturbations on Locally 
Rotationally Symmetric (LRS) class II spacetimes, exploiting a recently developed 
covariant and gauge-invariant perturbation formalism. From the Klein-Gordon 
equation and Maxwell's equations, respectively, we derive covariant and gauge- 
invariant wave equations for the perturbation variables and thereby find 
the generalised Regge- Wheeler equations for these LRS class II spacetime 
perturbations. As illustrative examples, the results are discussed in detail for 
the Schwarzschild and Vaidya spacetime, and briefly for some classes of dust 
Universes. 



PACS numbers: 04.20.-q, 04.40.-b 

1. Introduction 

The covariant 1+3 approach (see [1] and references therein) has proven to be a powerful 
tool in relativistic cosmology, especially through its application of the gauge-invariant, 
covariant perturbation formalism. This perturbation formalism works extremely well 
in cosmological applications when the background model is homogeneous and isotropic, 
that is of Friedman-Lemaitre- Robertson- Walker (FLRW) type. However, if the space- 
time under consideration has less symmetry, the 1+3 approach is no longer ideally 
adopted because its splitting in 'time' and 'space' relative to the fundamental observer 
is not sensitive to another preferred direction apart from 'time'. The description of 
space-time through covariant quantities, defined in the observer's restspace, is simply 
blind to a second distinguished direction. The 1+1+2 approach [2] remedies this 
by slicing the 'space' further into 'sheets' orthogonal to a second preferred direction. 
Analogously, the quantities of the restspace are further covariantly split in such a 
way that obtained quantities still have a clear geometrical or physical meaning. The 
1+1+2 approach thus naturally extends the 1+3 approach and keeps its benefits. 

In this article, we illustrate the 1+1+2 formalism with its application to scalar and 
electromagnetic perturbations on spacetimes which are locally rotationally symmetric 
(LRS) [3,4] possessing a continuous isotropy group at each point and hence a multiply- 
transitive isometry group. Since LRS spacetimes exhibit locally a preferred spatial 
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direction, the 1+1+2 formalism is ideally suited for a covariant description of these 
spacetimes, yielding a complete characterisation in terms of invariant scalar quantities 
which have either a clear physical or direct geometrical meaning. Such a covariant 
classification of LRS perfect fluid spacetime geometries has already been presented 
in [5], whereas orthonormal frame methods have been employed in [6-8]. We will 
include LRS imperfect fluids in our treatment but the emphasis will be mainly on 
LRS class II spacetimes, that is LRS spacetimes without vorticity terms such that the 
sheet becomes a genuine 2-surface. 

Schwarzschild black hole perturbations are well understood and it has been known 
for a long time that they are all governed by master equations known as the Regge- 
Wheeler equation [9], a Schrodingcr equation with a slightly different potential for 
scalar, electromagnetic and gravitational perturbations, respectively [9-14], due to 
the differing spins of the perturbing fields. Using the 1+1+2 formalism, we find the 
covariant generalisation of the Regge- Wheeler equation for scalar perturbations, as 
described by the Klein-Gordon equation, for all LRS spacetimes, and present the 
generalised Regge- Wheeler equation for electromagnetic perturbations, governed by 
Maxwell's equations, in the case of LRS class II spacetimes. We discuss the resulting 
wave equations in detail for Schwarzschild and Vaidya spacetimes [15, 16], the latter 
being closely related to the former by having a non-increasing mass. We also describe 
(source-free) electromagnetic perturbations on the Schwarzschild geometry by a linear 
first order system of ODEs plus an algebraic constraint, once spherical and time 
harmonics have been introduced. This allows for a quick determination of some 
electromagnetic field configurations, such as the solutions describing a static magnetic 
dipole or a static uniform magnetic field at infinity. We also discuss the key wave 
equations for two classes of dust Universe models. Specifically we have a look at 
the inhomogeneous Lcmaitre-Tolman-Bondi (LTB) spacetimes [15, 17-19], and the 
spatially homogeneous Kantowski-Sachs (KS) spacetime [15,20,21] in the spherical 
symmetric case. 

2. Preliminaries 

Let us briefly review the basics of the 1+1+2 covariant approach and introduce the 
required notation, referring the reader to [2] for a more elaborate exposition. We 
adopt units such that 8itG = 1 = c and Einstein's equations are G a b — T a b — kg a b- 

While the 1+3 approach [1] provides a threading of spacetime into 'time' and 
'space' with the means of a timelike unit vector field u a (u a u a = — 1) , representing 
the observers' 4-velocity, the 1+1+2 formalism in addition performs a slicing of the 
'space' with the help of a spacelike unit vector field n a , which is orthogonal to u a : 
n a n a — 1, u a n a = 0. The 1+3 projection tensor h b = g b + u a u b combined with n a 
gives rise to a new projection tensor N a b, 

N a b ee h b - n a n b = g b + u a u b - n a n b , (1) 

which projects vectors orthogonal to n a and u a (n a N a b = = u a N a b) onto 2- 
surfaces (N a a = 2) which are referred to as the 'sheets'. Its volume element, the 
alternating Levi-Civita 2-tensor, is derived from the volume element e a bc = u d rjdabc 
for the observers' restspaces by 

£ab = £ab c n c = u d rj dabc n c ; e ah n b = = e (a6 ) . (2) 
The covariant 1+3 threading irreducibly splits any 4-vector into a scalar part 
parallel to u a and a 3-vector part orthogonal to u a . Furthermore, any second rank 
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tensor is covariantly and irreducibly split into scalar, 3- vector and projected symmetric 
trace-free (PSTF) 3-tensor parts. The 1+1+2 slicing takes this split further by 
irreducibly decomposing the 3- vectors and PSTF 3-tensors with respect to n a . For 
example, any 3-vector tp a can now be irreducibly split into a scalar, 'J, which is the 
part of the vector parallel to n a , and a vector, ^5> a , lying in the sheet orthogonal to n a : 

^ a ^^n a + ^ a , where * = ip a n a , and V a = N ab ^ b = tp a , (3) 

where we use a bar over an index to denote projection with N ab . Similarly, any PSTF 
tensor, ip a b, can now be split into scalar, vector and tensor (which are PSTF with 
respect to n a ) parts: 

ipab = i\ab) = * (n a n b - \N ab ) + 2^ {a n b) + * afc , (4) 

where 





= n a n b i' ab = - 


-N ab ^ ab , 




= N a b n c ^ bc = 




*a6 


= (N {a c N b) d 


- \N ab N cd } 4> cd = * {ab} 



(5) 

We use curly brackets to denote the PSTF with respect to n a part of a tensor. Note 
that for 2nd-rank tensors in the 1+1+2 formalism 'PSTF' is precisely equivalent to 
' trans verse- traceless ' . 

In the 1+1+2 formalism three derivatives are used, which u a and n a define, for 
any object tp...'": 

^ a ... b c - d ^u e V^ a ... b c - d , (6) 
^ a ... b c - d =n e D e i.,-^^V'--V''/. C ---' l /Vi^, M , (7) 
5 e ^ a ... b c --- d = NjNj---N b 9 N h c ---N/D j i }f ... g h ^ . (8) 

The hat-derivative, analogously defined to the usual time (dot) derivative, is the 
derivative along the vector field n a in the surfaces orthogonal to u a . The 5-derivative, 
defined by equation © similarly to the spatial derivative D, is a projected derivative 
on the sheet, with projection on every free index. 

With these definitions the spatial derivative of n a orthogonal to u a is decomposed 
as follows: 

B a n b = n a a b + \4>N ab + ^e ab + ( ab , (9) 

where 

a a = n c B c n a = n a , (10) 
4> =S a n a , (11) 
C =\e ab 8 a n b , (12) 
Cab = S {a n b} . (13) 

These have a meaning analogous to the kinematical quantities 4-acceleration u a , 
expansion 9, vorticity io a and shear o~ ab of the congruence u a in the 1+3 formalism, 
which are related via 

V a Mf, = -u a ii b + \6h ab + e abc u c + a ab ; (14) 

more precisely, travelling along n a , <f> represents the sheet expansion, C ab is the shear 
of n a (distortion of the sheet), and a a its acceleration, while £ represents a 'twisting' 
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of the sheet (the rotation of n a [22]). The other derivative of n a is its change along 
u a , 

h a = Au a + a a where a a = h a and A — n a u a . (15) 

The new variables a a , 4>, £, Cab: -4 and a a are fundamental objects in the spacetime 
within the 1+1+2 approach since the spacetime geometry can be inferred from them. 
They are treated on the same footing as the kinematical variables of u a in the 1+3 
approach. 

The kinematical 1+3 quantities, the Weyl curvature and the energy- momentum 
tensor need also to be split. For the 4-acceleration, vorticity and shear, one arrives at 

u a =An a + A a , (16) 

uj a = Qn a + il a , (17) 

o- ab = S (n a n b - \N ab ) + 2E (a n b) + S a6 , (18) 
while for the electric and magnetic Weyl tensors one gets 

E ab = £ (n a n b - \N ab ) + 2£ {a n b ) + £ ab , (19) 

H ab = H (n a n b - ±N ab ) + 2H {a n b) + H ab ■ (20) 

The fluid variables heat flux q a and anisotropic pressure -K abl inferred from the energy- 
momentum tensor T ab together with the energy density \x and the isotropic pressure 
p, are split into 

q a = Qn a + Q a , (21) 
n ab = II (n a n b - \N ab ) + 2n (a n b ) + n ab . (22) 

Finally, the covariant 1+1+2 equations for all the above introduced quantities 
may be obtained from the Bianchi and Ricci identities for both the timelike unit 
vector field u a and the spacelike unit vector field n a . 



3. LRS class II spacetimes 

In this section, we discuss LRS class II spacetimes in terms of the 1+1+2 formalism. 
The discussion follows in parts van Elst & Ellis [5] but generalises their treatment 
of LRS class II perfect fluids towards imperfect fluids employing a new, somewhat 
streamlined notation. 

The 1+1+2 formalism is ideally suited to study spacetimes which exhibit local 
rotational symmetry (LRS) because they have a unique preferred spatial direction at 
each point, defined covariantly, for example, by an eigendirection of a degenerate rate 
of the shear tensor field or by a vorticity vector field. This direction constitutes a local 
axis of symmetry - all observations are identical under rotations about it and are the 
same in all spatial directions perpendicular to it [3,4]. Hence, after a 1+1+2 split, if 
the spacelike unit vector field n a (n a n a = 1, n a u a — 0) is chosen parallel to the axis 
of symmetry, all covariantly defined 1+1+2 vectors and tensors must vanish due to 
the LRS symmetry. It follows that LRS spacetimes may be covariantly characterised 
by scalar quantities (modulo equations of state for the matter variables), namely 

LRS: {A,e,(t>,^,E,Cl,S,H,fi,p,Q,Il,A} . (23) 

The variables l|23ll thus fully describe LRS spacetimes and are what is solved for in 
the 1+1+2 approach. For LRS class II, one requires in addition that the vorticity 
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terms vanish, Q = = £, which further constrains the magnetic Weyl curvature H. to 
vanish. Thus, a substantially smaller set of scalars describes LRS class II spacetimes: 

LRS class II : {.A, 6, 0, E, £, n,p, Q, II, A} . (24) 

These LRS class II quantities satisfy a bunch of covariant evolution and/or propagation 
equations, obtained from the Bianchi and Ricci identities for the unit vector fields u a 
and n a , respectively: 
Propagation: 

+ (!0 + S)(§0-£)-§( M - 
Q , 



s 

£■ 



¥ 

§A + §n = 

Evolution: 



f^E 



A)-s - |n 



(5 + in) + (±E - 10) Q ; 



(25) 
(26) 
(27) 



rfl - 



; n 



(e - §0) (a - \<j>) + q , 

i(/xH 

t0)n 



A^ + 2 



3p 



jE-0)£ + i(E- 



2A)-£ 
2 - 3 (/* H 



-in. 
p) (s 



Propagation/evolution: 
A- 6 = - (A + <f>) A + i6» 2 + §E 2 + i (fi + 3 P - 2A) 
/i + Q = + - (</>- 

Q + p + ti= -(§</> + A) n- 



2A)Q 

r6» + E) 



2 

|sn 



(28) 
(29) 
(30) 

(31) 
(32) 
(33) 



- {H + P)A 

This system of equations extends the corresponding system in [5] (see section 6) to 
imperfect fluids. 

Since the vorticity U, vanishes, the unit vector field u a is hypersurface-orthogonal 
to the spacelike 3-surfaces whose intrinsic curvature can be calculated from the Gauss 
equation for u a . For the intrinsic Ricci-curvature, one finds from this 

^2 1, 



3 R ab = [§ (M + A) + £ 

+ [| (m + a)- 



-i0E 



i£ 
2° 



ffl 2 ] ?i Q n b 
i^ + i^-f^] 



(34) 
(35) 



The intrinsic Ricci-scalar of the 3-surfaces is therefore 

3 R = 2 [fi + A - \6 2 + |E 2 ] . 

This relation constitutes the generalised Friedman equation. 

On the other hand, the additional vanishing of the sheet distortion, £, implies that 
the sheet is in this case a genuine 2-surface. The Gauss equation for n a together with 
the 3- Ricci identities determine the 3- Ricci curvature tensor of the spacelike 3-surfaces 
orthogonal to u a to be 



3 Rab = - [4>+\<t> 2 

thus giving for the 3-Ricci scalar 
3 R = -2 



n a n b 



K 



1 A? 



-K 



N ab 



(36) 



(37) 



Here, K is the Gaussian curvature of the sheet, 2 R a b — KN a b- Combining the last 
equation with equations l(35jl and 125(1 . we can express the Gaussian curvature K in 
the form 



K 



(p + A) — £ 



in + \<j? 



(38) 
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Evolution and propagation equations for K can be calculated from equations i(25|) - (|3"U)l 

to give 

K= - (§0-£)if , (39) 
K = -(j>K , (40) 
Every scalar t/j in LRS II has to satisfy the commutation relation 

^-ii = -Aijj + {\e + Y)i) . (41) 

We applied this commutator relation to the Gaussian curvature K in order to check for 
consistency of the equations ((25(1 - 133(1 . which is indeed guaranteed due to equations 
jSg> and gDJ). 

From equation 139|) follows the neat result that whenever the Gaussian curvature 
K of the sheet is constant in time but non- vanishing, that is, the sheets have either 
spherical or hyperbolic geometry, the shear E is proportional to the expansion 9: 

K ^ and K = => E = §0 . (42) 

As a matter of frame choice, it is always possible to pick an observer with k = 
0, corresponding to a 'static' observer. Adopting such a choice leads to great 
simplifications in the equations, the caveat being that such a choice is often not the 
most natural one, e.g., in Friedman-Lemaitre-Robertson- Walker (FLRW) spacetimes. 
For example, taking a static observer to describe spherically symmetric spacetimes, 
the constraint 1(38(1 allows for a covariant definition of the radial parameter r via 

r~ 2 = ! (/i + A) — £ — |ll + j(j) 2 ; r = = 6 a r . (43) 

(A non-static observer has r + 1 0). In situations with spherical symmetry it is 
convenient to express the hat-derivative in terms of the radial parameter r. If we 
associate an affine parameter p with the hat-derivative, we find using equation 140|) 

dX 1 dX 

x -~r^^ x ^ (44) 

where r denotes proper time and X may be any scalar. 

In summary, we have found that in the 1+1+2 formalism, all LRS II spacetimes 
are covariantly described by the scalar equations In this case, the sheets are 

genuine 2-surfaces whose extrinsic geometry is determined by the Gaussian curvature 
K, which can be computed from 1)38(1 . Furthermore, due to ((42(1 . we stress that 
static spacetimes whose sheet is either spherical or hyperbolic must have their shear E 
proportional to the expansion 9; more precisely, such imperfect fluid spacetimes fulfill 
the relation E — |# = 0, generalising the well-known fact that perfect fluid spacetimes 
within this subclass possess vanishing shear E and expansion 9. 



4. Perturbations on LRS class II backgrounds 

In this section, we investigate scalar and electromagnetic perturbations on LRS class II 
background spacetimes. The scalar end electromagnetic fields are treated as test fields, 
i.e., they do not affect the geometry of the background. Hence, the Stewart- Walker 
Lemma [23] ensures that these fields are gauge- invariant. Our goal is the derivation 
of the concomitant covariant and gauge-invariant wave equations which describe these 
perturbations. These wave equations are the generalisation of the well-known Regge- 
Wheeler equations from the Schwarzschild case towards LRS class II spacetimes, yet 
written in covariant notation. 
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4-1- Commutation relations 

In order to achieve this aim, we will make frequent use of commutation relations. All 
first order scalars ip have to satisfy the following commutation relations between the 



different derivatives introduced above: 

= - Aip + (§0 + £) i> , (45) 

5 J - N a b (8 b Tpy = -!(£- §0) 5 a iP , (46) 

6 a $ - N a b (S b TP) = + \<t>8 a iP , (47) 

S [a S b] i> = . (48) 

Analogously, the commutation relations for first order 2-vectors ip a read 

k-k = - Ma + {\0 + E) 4>- a , (49) 

Sa^Pb - N a c N b d (5 C M = - i (E - §0) S a ip b , (50) 

5 a ip b - N a c N b d (6 c ip d j = + i# Q V>6 , (51) 

5[ a 5 b ]ip c = - Kip[ a N b ] c . (52) 



4-2. Harmonics 

It is useful to expand all first-order perturbations in harmonics. Note that all functions 
and relations below are defined in the background only; we solely expand first-order 
variables, so zeroth-order equations are sufficient. 

In analogy to the spatial harmonics defined in [24], we introduce dimensionless 
sheet harmonic functions Q, defined on the background, as eigenfunctions of the 2- 
dimensional Laplace-Beltrami operator such that for positive, negative or vanishing 
curvature K 

k 2 

6 2 Q = — -Q , Q = = Q (0 < k 2 g R) . (53) 
The function r is, up to an irrelevant constant, covariantly defined by 

which, in the light of equations (|39|I - H40|) . is nothing but the covariant version of the 
common setting K = n/r 2 with k = ±1 for spherical and hyperbolic 2- geometry, 
respectively, where the constant is fixed by the relation However, the covariant 

definition l|54|l extends to the plane-symmetric case as well, where there is no natural 
length scale. We can now expand any first order scalar tp in terms of these functions 
formally as 

^ = E4 fc) Q (fc) =^sQ, (55) 

k 

where the sum over k is implicit in the last equality. The S subscript reminds us that 
tp is a scalar, and that a harmonic expansion has been made. 

We also need to expand vectors in harmonics. We therefore define the even 
(electric) parity vector harmonics as 

Qi fc ) = rS a Q^ => Q- a = Q = Qa, S 2 Q a = (l - k 2 ) r~ 2 Q a ; (56) 
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where the (k) superscript is implicit, and we define odd (magnetic) parity vector 
harmonics as 

Q<*> = re ab 5 b Q^ => 4 = = 4, b 2 Q a = (l - k 2 ) r - 2 Q a . (57) 

Note that Q a = s a bQ b Qa = ~£abQ b , so that e ab is a parity operator. The crucial 
difference between these two types of vector harmonics is that Q a is solenoidal, so 

S a Q a = , (58) 

while 

S a Q a = ~k 2 r- l Q . (59) 

Note also that 

e ab S a Q b = , and e ab 5 a Q b = +k 2 r~ 1 Q . (60) 

The harmonics are orthogonal: Q a Q a — (for each k), which implies that any first- 
order vector ip a can now be written 

= + ^l k) Q { a k) = ^Qa + Ma ■ (61) 

k 

Again, we implicitly assume a sum over k in the last equality, and the V subscript 
reminds us that iji a is a vector expanded in harmonics. 

We like to point out that the harmonics introduced here naturally generalise the 
spherical harmonics used in [2]. In particular, the various formulae for scalar and 
vector spherical harmonics stated in [2] also hold for our generalised harmonics. 

4-3. Scalar perturbations 

Let us consider perturbations of LRS spacetimes due to a massive scalar field, ip, with 
mass M, whose equation of motion is the Klein-Gordon equation, 

(ff Qb V a V 6 + M 2 ) j> = (V a V a + M 2 ) yj = . (62) 

We investigate this simple equation first in an arbitrary spacetime and specialising to 
LRS II spacetime afterwards. The corresponding 1+3 equation is easily derived using 
gab _ foab _ u a u b an( ^ following expression for the spatial Laplacian of a scalar 
field, 

D Q D Q V = h ab D a Dbi> = h ab V a V b Tjj + H , (63) 

and leads to the wave equation 

$ - D a D a V> + 9ip- " a DaV> + M 2 xjj = . (64) 

Thus, for a general spacetime, the evolution of a scalar field is affected by effects of 
expansion and acceleration, as measured by the fundamental observer. 

The 1+1+2 form of equation l|64|) is readily obtained using h ab — N ab + n a n b as 
well as the relation 

6 a 6 a ip = N ab 5 a 5 b T/j = 7V Qb D Q D b V - H> ■ (65) 
The resulting equation reads as 

^ - } + 9ip - (A + 4>)$ + (a a - A a - 6 a )S a ij} + M 2 V = . (66) 
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Equation (|66[1 is the fully split Klein-Gordon equation for a massive scalar field in 
an arbitrary spacetime, given in a covariant and gauge-invariant fashion. For LRS 
spacetimes, the wave equation <|66[1 closes, yielding 

$ - V> + 6ip - (A + <t>)i> + (M 2 - S 2 )ip = . (67) 

We emphasise that this is true for all LRS spacetimes. Equation H67JI constitutes the 
generalised Regge- Wheeler equation for scalar perturbations on all LRS background 
spacetimes. While the generalised Regge- Wheeler equation 1|67[) is not affected by 
rotational effects at all, it is in general sensitive to the observer's acceleration A, 
the spacetime expansion 8, the sheet expansion <j>, as well as the mass M of the 
perturbation field ip. 

For LRS class II spacetimes, it is sometimes convenient to rescale the scalar field 
ip according to ip = r -1 ^, where the function r is defined via l|54|) . In terms of the 
rescaled field VP, the Regge- Wheeler equation (|67|l reads 

§ - I + (E + §0) * - A^ + [M 2 - £ - \ [p - ip + 4A) - 5 2 ] * = . (68) 

This form of the Regge- Wheeler equation has several advantages over iffiTjl : firstly, it 
allows to introduce the Gaussian curvature K at the expense of the Weyl curvature 
£, for example, and secondly, it simplifies in vacuo. We emphasise that this form 
is generic: we will show in the next section that electromagnetic perturbations are 
described covariantly by the same equation but having a different potential [that is, 
the term in square brackets in (JHEJl] once appropriate harmonics have been used to get 
rid of the sheet Laplacian S 2 . 



4-4- Electromagnetic perturbations 

In accordance with ©, the electric, magnetic and current 3- vector fields are irreducibly 
decomposed into scalar and 2-vector parts as 

E a = Sn a + S a , B a — SSn a + 38 a , J a = J 'n a + J a . (69) 

Thus, Maxwell's equations for electromagnetic test fields, sourced by the total electric 
charge p c and the currents J* and J! a , respectively, on a LRS class II geometry become 
in the sheet approach (the splitting for arbitrary spacetimes is displayed in [25]): 

i + 5 a S a = - <\>S + p G , (70) 

,% + 8 a ,% a = -<f>&, (71) 

i - e ab 5 a ,98 b = + (E - §0) S _ J , (72) 

<% + e ab 5 a S b = +(E-§0)<?, (73) 

ia + e ab {$> - S b dj = - (\(j>. + A) e ab @ b - (|S + §0) S a - f a , (74) 

#a - s ab (i b - 5 b $) = + {\<t> + A) e ab S b - (|E + §0) <% a . (75) 

Note that e ab is a parity operator, s ac £ c b = —N ab , thus we can switch parity between 
2-vectors X a and Y a via 

X a = e ab Y b ^=>Y a = -e ab X b . (76) 

Using this relation, the parity-reversed form of equations ()74(l and l|75|l read 

i-a + s ab @ b -5 a S = - (\<t> + A)S a - (|S + §0) e ab m b , (77) 

@a ~ S ab i b - 6 a & = + (|0 + A) @ a ~ (|E + §0) Sab^ + E ab J b . (78) 
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From Maxwell's equations, one deduces (by a somewhat tedious calculation exploiting 
the earlier displayed commutation relations as well as the above parity-reversed 
equations) the following wave equations for the electromagnetic fields along the 
distinguished direction n a : 

i - i - (a + 2<f)) i - (s - §0) s 

- [S 2 + \<f -2£ + (±0 + E) (|E -B)-i(p-3p + 4A)] S 

= p c + J + {<t> + A) Pc + 0j? 1 (79) 

& - £ - (A + 24>) J - (s - f e) m 

- [5 2 + \4> 2 -2£+ (|0 + E) (|E - 0) - § (n - 3p + 4A)] St 

= . (80) 

Equations ((Til l!5U|) are the generalisation of the famous Regge- Wheeler equation, for 
electromagnetic perturbations on the Schwarzschild background, towards LRS class 
II spacetimes, although written in covariant manner. We emphasise that <§ and 33 
decouple from each other and that in the absence of sources the equations are identical 
closed wave equations. 

It will turn out advantageous to rescale the fields and the sources under 
consideration similarly as above for Klein-Gordon fields; that is, we define 

S = r- 2 E, 3§ = r- 2 B, Pc = r- 2 g c , J = r~ 2 J , (81) 

and substitute into the wave equations H79|) - (|80|l to obtain 

E - E - AE + (E + §0) E - [S 2 + (E - §0) (2E + |0)] E 

= g c + j + Ag c + (T, + l9)j, (82) 

B - B - AB + (E + ±0) B - [5 2 + (E - §0) (2E + \0)] B = Q. (83) 

Again, if we neglect the source terms, the equations become identical closed wave 
equations. Moreover, the equations H82(l - (|83|l are of the same form as equation (|68|l 
for the Klein-Gordon field, the only difference being a slightly altered potential term. 

The wave equations for the electromagnetic perturbations lying in the sheet are 
derived analogously and read 

4-4 - (3A + 0) i- a - (E - §0) 4 

- £ - A (A + (/>) + |0 2 - |0E - |E 2 + | {n - 3p + 4A) - S 2 ] S a 

2A) 8 a S + 3Ze ab <% b + ( 3^E -Q + O-A) e ah SS b 



- 5 aPc - J- a -{9- |E) f a , (84) 

S$- a -d- a - (3A + 0) J a - (E - §0) m- a 

+ [i0 2 - £ - A {A + (/>) + |0 2 - |0E - \Y? + i (fi - 3p + 4A) - S 2 ] SS a 

= {4> - 2A) b a SS - 3E£ Q6 <# b - (3^E -Q + §-Aj e ah S h 

+ e ab [5 b J ~ J b ~{ 1 1 cj ) + 2A) f h \ . (85) 

In contrast to equations (|79() - (|80|l , these do not decouple-not even in the absence of 
sources. For example, in addition to the source terms, the magnetic 2-vector field 3S a 
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gives rise to forcing terms for the electric 2-vector field & a , which is also forced by the 
'radial' electric field. 

However, it turns out that in the absence of sources knowledge of the 'radial' 
part $ (or BS) of the perturbations, e.g., by solving the concomitant Regge- Wheeler 
equation l|82[l . suffices to completely determine the electromagnetic perturbations. To 
see this, we expand all perturbations into harmonics and decompose the governing 
Maxwell's equations into their harmonic components. (For details, we refer the reader 
to Ref. [2], Section E). From the scalar equations (|70|I - (|73JI we find for each fixed 



harmonic index k 

is + <Ks = + , (86) 

,#s + <^s = + , (87) 

A - (E - §0) Ss = + , (88) 

#s - (E - f 0) #s = - ■ ( 89 ) 



It is thus obvious that a solution for the 'radial' fields S and SB determines the sheet 
fields S a and 38 a , respectively. 

The 2-vector equations (|74|) and (|75|) . or l)77|) and J7SJ), respectively, split into 



two parts; the first part is the odd parity part, given by 

'iy + - = - {\4> + A) SSy - (iE + §0) Jv , (90) 

&\i-£m + \£s = + {±4> + A) S-y - + 16) J v , (91) 

and the second part is the even parity part, given by 

i v - ®m = +(\<l> + A)Mv- (fE+§0)<?v , (92) 

Jv + ^v= - (\<i) + A)i\i- (iE+|6»)^v (93) 



We remark that the even parity equations are redundant since propagating the 
constraints l|88fl - H89|) and inserting them into l|92|) and (|93|l just yields back 
equation (|86|) and l|87l) . On the other hand, the odd parity equations are not implied 
by the scalar ones and have been used in the derivation of the generalised Regge- 
Whceler equation l|79|l . If one eliminates i\j and 3§y from equations Ij92(l - I|93|l using 
the constraints (|88|I - H89|I . it becomes obvious that the electromagnetic perturbations 
fall into two distinct classes whose equations decouple from each other, namely: 

polar perturbations : {Ss, £{/, , (94) 

axial perturbations : {^S)^v,^v} ■ (95) 

Moreover, the resulting equations involving either Ss and Sy or 3§s and are 
identical. 

5. Examples 

We have shown that scalar and vector perturbations in LRS class II spacetimes are 
governed by simple master equations. We shall now discuss these master equations 
in some specific spacetimes. First, we shall recover the well known Regge- Wheeler 
equations for the Schwarzschild spacetime, which we then generalise to the Vaidya 
radiation spacetime. We then discuss solutions from the Lemaitre-Tolman-Bondi and 
Kantowski-Sachs families. 
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5.1. The Schwarzs child spacetime 

5.1.1. The background Schwarzschild spacetime is fully determined by any two of 
three non-zero scalar functions <j> , A , £. These functions obey the background 
equations 

$ = - \<t> 2 + M , (96) 
A= -A{<t> + A) ; (97) 

together with the constraint 

e + A(j) = . (98) 

The parametric solutions for these variables are [2] 

e = ~, (99) 



2 / 2m , 

= -Jl , (100) 

r V r 

/ \ -1/2 



r 



A = -s 11- — ) . (101) 



These form a one-parameter family of solutions, parameterised by the constant m, 
which is just the Schwarzschild mass. The (exterior) Schwarzschild solution is given 
for 2m < r < oo. 



5.1.2. The perturbations In order to make the connection with the standard Regge- 
Whceler equations [13], it is instructive to contrast the found covariant and gauge- 
invariant wave equations for scalar and electromagnetic perturbations with the ones 
governing gravitational perturbations. The fundamental object for the study of 
gravitational perturbations on a Schwarzschild background is the Regge- Wheeler 
tensor W a b [2] , defined by 

W ab = i0r 2 Cafc - ir 2 £-M { A}£ , (102) 

a dimensionless, gauge-invariant, transverse-traceless tensor which obeys the wave 
equation 

W {ab} - W {ab} - AW {ab} + 4> 2 W ab - S 2 W ab = . (103) 

If one expands equation i|103|) m f° spherical tensor harmonics, the odd and even 
parity parts of equation (|103fl both become 

L 



Wt-Wt- AW t 



, , 3£ 



W T = , (104) 



where Wj = Wj are the tensor harmonic components of W ab , £ — 1,2,... and 
L = 1(1 + 1) = k 2 . It turns out that equation p()4|l is actually the Regge- Wheeler 
equation [9] when written in appropriate coordinates. Converting from p, the affine 
parameter associated with the hat-derivative, to the parameter r, p — > r, and then to 
the 'tortoise' coordinate of Regge and Wheeler, 



2mln(--l) , (105) 
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and also introducing the Schwarzschild time via dr = \/ 1 — =2* dt, we find that l|lU4fl 
becomes: 



L 6m 



(106) 



where ~fr is the Regge- Wheeler potential for gravitational perturbations. 

It is now straightforward to obtain the familiar Regge- Wheeler equations for 
the case of scalar and electromagnetic perturbations. From equation (|68|l one finds 
immediately for a massive scalar field ijj with mass M, ip = r -1 ^/, that 

* - * - - [£ - M 2 + 5 2 ] * = . (107) 
Introducing scalar spherical harmonics and performing as above one readily gets 



d 2 d 2 , \ , / 2m 



dt 2 dr 



*s = , f s =[l- 



r 



L 2m o 
— + — + M 



(108) 



where is a scalar harmonic component of and ~fs is the Regge- Wheeler potential 
for scalar perturbations. This equation was originally derived in [10] employing 
the Newman-Penrose formalism. Notice the striking similarity between the Regge- 
Wheeler potential for the massless scalar field perturbations, and the potential for 
gravitational perturbations, ~fr [13]. 

The electromagnetic case is even simpler. For this case, the covariant Regge- 
Wheeler equation l|82|l for the electric field perturbation S ', S — r~ 2 E, reduces in the 
Schwarzschild background to 

E — E — AE - 8 2 E = . (109) 
Adopting once more scalar spherical harmonics and proceeding as before leads to 



d d ,„ \ „ „ ,, I _ 2m 



(110) 



where Yy is the Regge- Wheeler potential for electromagnetic perturbations (the 
magnetic case is completely analogous). The Regge- Wheeler equation for radial 
electric perturbations (|110f> . and its magnetic counterpart, was originally derived 
in [11] employing the Newman-Penrose formalism. 

It is remarkable that the (source-free) electromagnetic perturbations can also be 
described by a linear first order system of ODE's after expanding Maxwell's equations 
into spherical and time harmonics [2] (possible because the background is static). The 
time derivatives of first order quantities are decomposed into their Fourier components 
by assuming an e luJT time dependence for the first order variables; factors of iui just 
represent time derivatives, d/dr. Note that 

2m \ 2a 



uj = -Auj u = a 1 =— , (111) 

\ f J (pr 

arising from the commutation relation between the dot- and hat-derivatives. The 
harmonic function uj is defined with respect to proper time r of observers moving 
along u a , while a is the constant harmonic index associated with time t of observers 
at infinity. They are related by lot — at. Inserting time harmonics into equations 
lISBl l|HT|) one finds that equations (JHSJ-JS^J turn into constraints for the odd modes 
S\j and ^v, respectively, 

- iuor - iujr „ 

v = T s ' v = ~~ s ' (112) 
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which may subsequently be eliminated from equations (|90|I - H91|) . The final result is a 
linear first order system of ODE's, 

Y = AY , (113) 

where the matrix A is given as 

A -(i-~** 4 -a) 

and the perturbation vector Y is defined in the case of polar and axial perturbations 

by 

''polar = ( J? J and Y axia i = ( ^ S ) , (115) 



respectively. We emphasise the observation that in a stationary situation, e.g. u> = 0, 
the constraints (|112|) imply Sy — SSy — in the absence of charges, which means that 
electromagnetic perturbations cannot have a solenoidal part in the sheet. A further 
implication is that a pure magnetic field can only exist on the Schwarzschild geometry 
if it has no solenoidal contribution in the sheet and is stationary and/or sheet-like. In 
the stationary case, the system l|113|) can be solved analytically. For example, picking 
L = 2 in l|114f> easily yields the well-known solutions of a static magnetic dipole or a 
static uniform magnetic field at infinity (compare also [26-28]). 

5.2. Vaidya Spacetime 

5.2.1. The background The energy-momentum tensor of Vaidya's radiating sphere 
spacetime is that of a radiation fluid, e.g., 

T ab = tik a k b , k a k a = . (116) 

A convenient form of its metric for the case of outgoing radiation is given in terms of 
retarded time u = t — r* by 

g ab = -(l- ^M) du 2 - 2dudr + rW , (117) 



where m(u) is a arbitrary non-increasing function, m'{u) = dm(u)/du < 0; if m(u) 
is constant, the metric is equivalent to the Schwarzschild metric. Thus, the Vaidya 
spacetime describes a collapsing star that irradiates a part of its mass, where the total 
power output measured at infinity is given by -m'(tt) = 4?rr 2 /x [29-34]. 

A covariant description of the Vaidya spacetime is obtained as follows. Spherical 
symmetry requires the null- vector k a to be of the form k a = u a ± n a , where a plus 
(minus) corresponds to outflowing (inflowing) radiation. Thus, a static observer will 
encounter the radiation fluid as an imperfect one, with the energy-density fi, the 
isotropic pressure p = |/^, the radial heat-flux Q = ±/i and the trace- free part of the 
anisotropic sheet pressure II = Working in the static frame (r = 0), the non- 
stationary spherically symmetric Vaidya spacetime in the case of outgoing radiation 
is described covariantly by the following set of equations: 

A- 9= -A{A + <])) + 6{0-<i>) , (118) 

4> = -cj){\(t>-A-2e) , (119) 

<j> = -0(f> , (120) 
6 + 9 = -9(±(f> + 39 + 3A) . (121) 
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The corresponding constraints, implied by our frame choice , are 

K = r- 2 = ^ 2 ~£, (122) 

£ = §0 , (123) 

H = - 6(f) = 3p = Q = §n , (124) 

£ = fx-A(f) . (125) 

Note that the equations governing the Vaidya spacetime reduce to the ones describing 
Schwarzschild spacetime [see (|96[) -(l98 fl ] when the expansion 9 vanishes. 

In order to solve the equations, we find it useful to express our derivatives in terms 
of coordinates {r, u}, where r is the radial parameter and u labels time. In particular, 
we have p — p(r, u) and t = r(r, u) and the derivatives become for an arbitrary scalar 
X: 

r6dX dX , x 

Fix 

<I27> 

One neat possibility to define the coordinate u is through the relations 

1 2 

u = -it = - = — , (128) 
r r<p 

which is compatible with the commutator relation l|45|l and identifies u with the it- 
coordinate used in the metric 1117JI . as will become clear in an instant. Upon switching 
to the coordinates {r, u} and putting A = Ascb — & [as suggested by equations (|1 lfJfl 
and l|121(l ]. the system of equations transforms into 

(\r 2 4> 2 d r - 8 U ) Ach = - \r<t> [Ascu (Ascu +<b + 0) + \<j>0] , (129) 
rd r (j> = -0 + 2Ach , (130) 

2d u cb = -r9<j) 2 , (131) 

r<t)d r 6 = - ((/> + 6Asch) ■ (132) 

The solutions may be readily found by making an educated guess. Clearly, by spherical 
symmetry, the radiation has to move outward radially and the expansion <j) of the 
radial congruence n a therefore should be of the same form as for the Schwarzschild 
spacetime but with a decreasing mass parameter. Since equation (|130fl is identical 
with the analogous one for the Schwarzschild spacetime, we are led to identify -4sch 
with the acceleration of a static Schwarzschild observer, in particular, ,4sch represents 
the acceleration caused by the star's instantaneous gravitational field. It is then 
straightforward to work out all other quantities. The result is: 



r 



e = U = §s t (134) 

A = Ac - 0; As* = ^ (l - ^) , (135) 

2mw 1 _^r 1 =3p=Q=|II) (136) 



2m{u) 



(137) 
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It is worth pointing out that these solutions can also be worked out directly from the 
Vaidya metric l|117fl by choosing the splitting unit vectors to be 

2m(u)Y 1/2 ( d 

and 



"* = I 1 - I I g-J (138) 



du J \ r J \dr 

respectively. 

The obtained solutions may be interpreted as follows. First note that from 
equation i|13ti|) follows that a physically acceptable, that is, positive energy density 
\i requires a non-increasing (but otherwise unrestricted) mass, m!{u) < 0, suggesting 
that a part of the star's mass is released as radiation during collapse, with total 
luminosity — m'(u) = Anr 2 ^, at infinity (in proper units). Further, the expansion 
9 and shear E are negative since m'(u) < 0, meaning that the spatial sections 
of the spacetime are contracting, which tends to push nearby observers together. 
Equation l|135|) hence tells us that an observer has to balance its always radially 
outward directed acceleration according to the diminishing gravitational attraction of 
the collapsing star and the growing contraction in order to stay at rest. Finally, the 
expressions for sheet expansion <fi and tidal forces £ are as for a star in equilibrium 
but with a time-dependent mass due to mass-radiation conversion during collapse. 



5.2.2. The perturbations With the background solutions for the Vaidya spacetime 
at hand, it is now an easy task to rewrite the covariant Regge- Wheeler equations 
for the case of scalar and electromagnetic perturbations into in terms of coordinates. 
From equation (|68|l one obtains immediately for a massive scalar field -0 with mass 
M, ip = r^ 1 *, that 

* - * - + 0* - [£ - M 2 + S 2 ] * = , (140) 
while equation l|82l) for the electromagnetic perturbations, S — r~ 2 E, gives 

E-E-AE + 6E-5 2 E = 0, (141) 

Introducing scalar spherical harmonics and adopting the coordinates {r, u} from above, 
we find the corresponding Regge- Wheeler equation 



2m(u) 



dr 2 



d 2 2m(u) d 



dudr 



dr 



L , 2m(u) 
+ A ^ 



AM' 



= 0, (142) 



where 3^s denotes the scalar harmonic component of either ^> or E and L = £(£ + 1); 
the parameter A takes the value one for perturbations in a scalar field and the value 
zero for electromagnetic perturbations. 

It is instructive to compare the equations l|142f> with the corresponding equations 
H1Q8(I and l|110(l in the Schwarzschild case. We found it most favourable to do the 
comparison by employing coordinates {t,r^} which are defined for r > 2m(u) by 



t = u + r* 



r* = r + 2m{u) 



In 



2m(u) 



(143) 
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In terms of the coordinates {i,r*}, the equations (|142H now become 

4m' (u) 



d 2 d 2 
{-l-h{u,r)] — + [l-h{u,r)} — , 

ot A or% r — 2m{u) \ot or. 



1 



2m(u) 



L , 2m(u) , , r9 
^7 + A — ^ + AM 



where the function h(u,r), 

h(u, r) = Am! (u) 



In 



2m(u) 



-1-1- 



^ s = 



2m(u) 



(144) 



(145) 



was introduced for brevity's sake. Clearly, when the mass m(u) stays constant, the 
equations (1 14411 reduce to the familiar Regge- Wheeler equations I|1U8|) and 1)1 10|) . 
respectively. Note that all new terms in the Regge- Wheeler equations for the Vaidya 
spacetime 111440 are proportional to the mass change in time, m'(u), and thus 
proportional to the expansion 9 [cf equation (|134|l ]. This was to be expected as, 
contrasting the Schwarzschild case, the only new ingredient in the covariant Regge- 
Wheeler equations for the Vaidya spacetime [see (|140|) and (|14H ] is proportional to 
the expansion 8 as well. Finally, the accompanying Regge- Wheeler potentials retain 
the Schwarzschild form but become time-dependent in the Vaidya case. 



5.3. More general spacetimes 

It was shown in [4] that the metric of every LRS class II spacetime can be given in 
diagonal form if local comoving coordinates are chosen: 

ds 2 = -A- 2 (t,x)dt 2 +B 2 {t,x)dx 2 +C 2 (t,x) [dy 2 + D 2 {y,k)dz 2 ] , (146) 

where D(y, k) = (smy , y , sinh y) for k = (1 ,0 , —1) labelling the closed, flat or open 
geometry of the sheet. The kinematical quantities are now 



(147) 
(148) 
(149) 
(150) 



4> = 


2 £ 

C ' 


A = 


A 

~ A ' 


6 = 


§+' 


E = 


2 ( B 

3 I B 



and equations l|39|l and l|4(JII can therefore be integrated to yield the Gaussian curvature 

K -c$T)' (151) 

where the integration constant C\ may be normalised such that C\ — k. From equation 
H2(jf> or (|28|) one immediately obtains for the heatflux 



(152) 
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while equation l|31|) directly yields the relation 
§ (/i + 3p - 2A) = - 

which can be employed to constrain the metric functions such that, e.g., /i + Zp > 
holds. If a simple equation of state, p — afi, say, is assumed, one finds from this 
relation an expression for the energy density fx. On the other hand, the field equations 
imply 



\i + A - K = 




(154) 



which is in general different from the expression for /i obtained from equation l|153fl 
and therefore gives an additional constraint for the metric functions. Finally, the 
pressure variables and the electric Weyl curvature can be calculated similarly but the 
resulting expressions are somewhat long and will not be stated here. 



A n A\ B C AC 

I" 2 I + b +2 c + 2 ac < (153) 



5.3.1. The perturbations When local comoving coordinates are chosen, the wave 
equations governing the scalar and electromagnetic perturbations in the sourcefree 
case are best given in form of equation (|67|l and (|82|l . respectively, since there only 
the kinematical variables enter and these are simple expressions in terms of the metric 
functions [cf. equations l|147[) - i|15UI) ]. The corresponding wave equations read in the 
case of scalar field perturbations ip as 







(155) 



while in the case of (rescaled) electromagnetic perturbations E (or B) they become 

2 




3^-2 (£) 
BC \C 



E = 0. (156) 



Once the metric functions A, B and C in (|146|l are given for a known LRS class 
II solution, these wave equations may readily be transformed into their concomitant 
coordinate analogues. 

Observe that the 'radial' parameter r, covariantly defined by the relations l|54|) . 
agrees with the metric function C(t, x) in the line element l|14tj|) in account of l|147fl 
and l|149[l - (|150f) . Thus the physical electromagnetic perturbations § are obtained 
from the rescaled ones E by the transformation <§ — C~ 2 (t, x)E. If the scalar field ip 
is analogously rescaled as -0 = C -1 (i, x)^, then the wave equation for the scalar field 
perturbation (|155fl takes the Regge- Wheeler form 



2 C C BC AC 2 
+ C C BC AC 



* = 



(157) 



It is thus obvious that once a harmonic decomposition has been applied to the 
equations l|156l) and (|157fl . the resulting wave equations will only differ in their 
corresponding potential term. 
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If in addition to the LRS symmetry further symmetries are present, then the wave 
equations for the perturbations simplify considerably. For example, for stationary LRS 
class II spacetimes we have 



* - <f + I — I * 



E-E 




M l + — 



C AC 



C AC 



* = 



E - 5 2 E = 



while for spatially homogeneous spacetimes (where one may set A = 1) we have 



*-*+§* 



E-i + | A. 



_ ,o C BC _n 

M -C-BC' 5 



* = , 



3^-2 (£) 
BC \C 



E = 



(158) 
(159) 

(160) 
(161) 



5.3.2. Application - LTB dust Universe As an example of an inhomogeneous 
spacetime, let us consider the Lemaitre-Tolman-Bondi (LTB) [15,17-19] dust Universe 
(with A = 0), whose metric is given by 

<Y') 2 

K 1 dr 2 + Y 2 [d-d 2 + D 2 (i),k)d<p 2 } . (162) 



-dr 



k-sf 2 (r) 

Here, a prime means d/dr and a dot will denote d/dt. Moreover, Y = Y(t,r) is 
the solution of Y 2 — 2m(r)/Y = — e/ 2 (r), wherein e = (—1, 0, 1) corresponds to the 
hyperbolic, parabolic and elliptic solution, respectively. The constraint equation can 
be integrated completely, which yields an additional function is(r). Therefore, there 
are three functions which can be prescribed at will: the 'mass' m(r), the 'energy' f(r) 
and the 'bang time' ts{r). 

It follows from equations (|147[l - l|150|l that the acceleration A has to vanish and 
that the non-zero dynamical variables take the form 



Vk^Tf 2 




(163) 
(164) 

(165) 



An expression for the energy density /i can be gained from equation 1)153(1 and yields, 
using the above mentioned constraint, 

2 m' , 

Finally, the wave equations (|156fl and (|157fl governing the perturbations 3? can be 
written as 



k-sf 2 
(Y'f 



Y" eff 
Y> k- ef 2 



^ 0> + {V -6 2 )0> = ; (167) 
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the potential V is given by 

YY' (y\ 2 

V = V EM = 3 — -2\-j (168) 
in the case of electromagnetic perturbations, whereas 

denotes the potential in the case of scalar field perturbations. 

5.3.3. Application - Kantowski- Sachs dust Universe As a further example we 
consider a spherically symmetric Kantowski-Sachs dust Universe with A = [15,20,21] . 
The metric writes as 

ds 2 = -dt 2 + B 2 (t) dr 2 + C 2 {t) dn 2 , (170) 

where the metric functions are conveniently expressed in terms of a parameter rj(t) 
satisfying dt — 2Cdr], namely 

B = m (77 tan 77 + 1) + 6 tan 77 , (171) 
C = c cos 2 77 ; m ,b , c = const . (172) 

Since this spacetime is spatially homogeneous, the acceleration A and expansion (f> 
of the fundamental observer's congruence have to vanish [cf. equations (|148|l l|147Jl ]. 
The remaining quantities are found to be 

, ~ , (173) 
777 cos 2 77 -3£ 

- ZBC 2 cos 2 n ' [ ' 

m (ri — 3 sin 77 cos 77 — 4 77 sin 2 77) + 6(1 — 4 sin 2 77) 
= c < 2 ^ c2 ^ * 11 , (175) 

m (77 + 3 sin 77 cos 77 + 2 77 sin 2 77) + b (l + 2 sin 2 77) 
£ = c — ii ! ^ 2 '1 V '1 . (176) 

The equations for the perturbations & now become in terms of the parameter 77 

id 2 AC 2 d 2 B (1 + 2 sin 2 77) -m cos 2 77 d , , 1 , 

I a?7 z ij z c*r^ is sin 77 cos 77 077 I 

where 

9 2_B — to cos 2 77 

V = V S = M 2 + — — g-i 178 

2BC 2 - cos z 77 

in the case of scalar perturbations, and 

3 m cos 2 77 — B (3 + 4 sin 2 77) 

V = = ' 2 V 2 ^ 179 

2BC 2 cos 77 

in the case of electromagnetic perturbations, respectively. 
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6. Conclusion 

Employing the 1+1+2 formalism of Clarkson & Barrett [2], we presented a covariant 
description of LRS class II spacetimes in terms of scalar quantities, which all 
have either a clear physical or geometrical meaning. We investigated scalar and 
electromagnetic perturbations (test fields) on LRS class II spacetimes and found that 
they are governed by covariant wave equations [see equations ijrJSj) and (JH2J], which 
are the covariant generalisations of the Regge- Wheeler equation, known to describe 
perturbations of the Schwarzschild spacetime. In particular, it was shown that both 
scalar and electromagnetic perturbations (in the absence of sources) are governed by 
master equations of the same form, namely the covariant generalised Regge- Wheeler 
equation, 

& - & - A & + (S + §0) & + (V - S 2 ) & = , (180) 

where the potential V is different for the two cases considered. To arrive at this 
specific form, one had to rescale the perturbations with the 'radial' parameter r, which 
is induced from the Gaussian curvature K of the sheet and can thus only be defined 
covariantly for LRS class II spacetimes [cf (|54() ]. The findings have been discussed in 
some detail for the particular cases of Schwarzschild and Vaidya spacetimes, and for 
some dust Universe models. While the master equations may be written simply as 
equation l|180|l in covariant form, they can become very untidy when written explicitly 
in coordinates, demonstrating some of the advantages of using a covariant approach. 
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